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ABSTRACT 

The inverse Compton scattering of interstellar photons off cosmic-ray electrons seems to play a more important 
role in the generation of diffuse emission from the Galaxy than thought before. The background radiation field 
of the Galaxy is highly anisotropic since it is dominated by the radiation from the Galactic plane. An observer 
in the Galactic plane thus sees mostly head-on scatterings even if the distribution of the cosmic-ray electrons is 
isotropic. This is especially evident when considering inverse Compton scattering by electrons in the halo, i.e. the 
diffuse emission at high Galactic latitudes. 

We derive formulas for this process and show that the anisotropy of the interstellar radiation field has a sig- 
nificant effect on the intensity and angular distribution of the Galactic diffuse 7-rays, which can increase the 
high-latitude Galactic 7-ray flux up to 40%. This effect should be taken into account when calculating the Galac- 
tic emission for extragalactic background estimates. 

Subject headings: cosmic rays — diffusion — Galaxy: general — ISM: general — gamma rays: observations — 
gamma rays: theory 



1. INTRODUCTION 



Recent studies show that inverse Compton scattering (ICS) of interstellar photons off cosmic -ray (CR) electrons seems to play a 



more important role in generat ion of diffuse 7-ray emission from the Galaxy than thought before (e.g., Porter & Protheroe 1997 
Pohl & Esposito 1998| , |VIS98b ). The spectrum of Galactic 7-rays as measured by EGRET shows enhanced emission above 1 GeV 



in comparison with calculations based on locally measured proton and electron spe ctra assuming the same spectral shape over the 
whole Galaxy ( |Hunteret al. 1997| , |Gralewicz et al. 1997| , |Mori 199% [SM97| , |MS98a| ). The 7-ray observations therefore indicate that 
their spectra on the large scale in the Galaxy could be different. Harder CR spectra could provide better agreement, but the 7-ray data 
alone cannot discriminate between the 7r°-decay and inverse Compton explanations. The ha rd electron spectrum h ypothesis seems 
to be more likely due to the probably clumpy distribution of electrons at high energies (e.g., ^ohl & Esposito 1998 ), while the hard 



nucleon spectrum can be excluded at the few a level on the basis of antiproton data above 3 GeV and positron measurements ( |MS98b 



MSR98). If the interstellar electron spectrum is indeed harder than that measured directly, then the ICS is a main contributor to the 



diffuse Galactic emission above several MeV and can account for the famous 'GeV excess'. 

The origin and spectrum of the Galactic diffuse 7-ray emission has a strong impact also on the extragalactic studies and cosmo- 
logical implications. The origin of the trul y extragalactic 7-ray ba ckground is still unknown. The models discussed range from 
the primordial blac k hole evaporation (e.g., Page & Hawking 1976 ) and annihilation of exotic pa rticles in the early Uni verse (e.g., 
Cline & Gao 1992) to the contribution of unresolved discrete sources such as active galaxies (e.g., Sreekumar et al. 1998), while the 
spectrum of the extragalactic emission itself is the main uncertain component of these models. The latter can be addressed only by 
the accurate study of the Galactic diffuse emission. Additionally, it is relevant to dark matter studies such as the search for signatures 
of WIMP annihilation, which is also a potential source of 7-rays, positrons, and antiprotons in the Galactic halo (see J ungman, 
Kamionkowski, & Griest 1996 for a review). 

The calculations of the ICS in most studies are based on the formula obtained by Jones ( 1 968) which remains a good approximation 
when considering isotropic scattering. However, the interstellar radiation field (ISRF) (apart from the cosmic microwave background) 
is essentially anisotropic. The stellar and dust emission originates in the Galactic plane and preferentially scatters back to the observer 
also in the plane. Therefore, the observer sees mostly head on collisions even if the e lectron distribution is isotropic. Though 
qualitatively such an effect was pointed out many years ago (e.g., Worrall & Strong 1977), the magnitude of this effect has not been 
evaluated accurately in the literature so far. Since ICS is one of the most important components of the diffuse Galactic 7-ray emission 
its accurate calculation is of interest for astrophysical studies. 

We have developed a propagation code ('GALPROP') which aims to reproduce self-consistently observational data of many kinds 
related to CR origin and propagatio n: direct measurements of nuclei, antiprotons, electrons and positrons, 7-rays, and synchrotron 
radiation ( |M9§ |MS98a| , |MSR98[ ). In this paper we evaluate the effect of anisotropy in the ICS of CR electrons off interstellar 
photons and present c orrespon ding formulas, which serve a basis for our detailed analysis of the Galactic diffuse emission in an 
accompanying paper ( SMR99 ). We make calculations with a realistic ISRF and also a simplified distribution of photons (thin 
emitting disk), which allows for simpler analytical treatment. The latter shows the effect without the uncertainty introduced by 
calculation of the real ISRF (which is itself a complicated task). For interested users our model including software and result datasets 
is available in the public domain on the World Wide Web[]. 

Throughout the paper the units h = c = m e = 1 are used. 

'http://www.gamma.mpe-garching.mpg.de/~aws/aws.html 
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ANISOTROPIC INVERSE COMPTON SCATTERING IN THE GALAXY 



2. INVERSE COMPTON SCATTERING IN ANISOTROPIC PHOTON FIELD 



We start from the general formula for the rate of photon-particle interactions ( [Weaver 1976| ): 

R = n e n 1 J dp-, J dp e / e (/ e )/ 7 G? 7 )^^a-(/> 7 ), 



(1) 



where « e ,« 7 are the electron and photon number densities, respectively, p e ,p-, are the momenta, / e (p e ),/ 7 (p 7 ) are the corresponding 
distribution functions in the laboratory system (LS) (the normalization being chosen so that J f en (p e .-,)dp e .~, = 1), 7 is the electron 
Lorentz factor, a is the cross section, and prime marks the electron-rest-system (ERS) variables. In the following treatment the 
electron distribution is assumed isotropic. 

2.1. Scattering off a single electron 

For an electron energetic enough the incoming photons in the ERS are seen as a narrow beam, ~ I/7 wide, in the forward direction. 
Therefore, we adopt Jones' (1968) approximation considering the incoming photons as a unidirectional beam in the ERS. Relaxing 
this assumption leads to overcomplexity of the final formulas and their derivation becomes lengthy and laborious while has little 
effect on the final result. An equivalent approach is valid for the LS, the photons are most ly scattered forward in a cone of ~ I/7 
wide and can be also considered as unidirectional. The latter will be used in Section 2.2.1 to set up kinematic limits on the angles 
involved. 

First consider monoenergetic electrons with distribution 



fe(Pe) : 



1 



;5{p e -p), 



A-Kp 2 e 

where 8{x) is the Dirac delta function. In terms of the electron Lorentz factor this gives (f2 e is the solid angle) 

1 



47T7 



2 5(71-7), 



where p = (3j, and f3 = 1 is the electron speed. Consider also monoenergetic target photons with distribution given by 

/ 7 (e~,n~) = 2 7 (0»4-^ e 7~ £ i)' 



(2) 



(3) 



(4) 



where e 7 = p~, is the photon energy, Q-,(9, </>) is the angular distribution of photons ( J dfl~, (2 7 = 1) at a particular spatial point, and 
(9, 4>) are the photon polar angle and the azimut hal angle, respectively ( Fig, nil. 
The Klein-Nishina cross section is given by (Jauch & Rohrlich 1976) 



do 



de' 2 d cost]' 



-e'[(l - cos rf) 



(5) 



where r e is the classical electron radius, e\ , e' 2 are the ERS energies of the incoming and scattered photons, respectively, and rj' is the 
scattering angle in the ERS. 

The upscattered photon distribution over the LS energy, £2, as obtained from eq. ([[]) is 



dR 

d€2 



J d cos f]' J deydCly J d^ ] dQ e f e (j ] ,Q e )f 1 (e 1 ,n j )e 2 1 ^-^j(^j — 



da 



dcosrj' 



(6) 



where J(e 2 /e2) = ^'2/^2 is the Jacobian (ERS — > LS), e' 2 is a function of cos 77' and LS variables, p' 7 = e^, and we have put n c . 7 = 1. 
Taking into account the 5-functions in eqs. (||),(Q), the integration over afe 7 ,£/7i,flff2 c , is trivial, and yields 



dR(^,e{) irrl 



de 2 



7 £ 1 £ 2 



c/cos??' / cm 7 2 7 (fi 7 )^- 4 + -r- sin V U 4"-; 7 



;'j(1 -cos 77') 



(7) 



where e' 2 = £2/[7(l -/3cos?/)]> and ?/ = it for forward scattering. 
The integration over d cos ?/ is also trivial (7 > £2) 



d€2 



ei(7-£2) 2 



dn 7 g 7 (^ 7 ) 



2-2 



£2 / 1 



where 



e 2 <27e , i/(l+2e' 1 ), e[ = ei 7 (l +/3cosC), 



(8) 



(9) 
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FIG. 1. — Angles involved in the ICS process: (8,ip), polar and azimuthal 
angles of the incoming photon, (K,ip), the same for the electron, R s , observer's 
position. See text for more details. 




FIG. 2. — The spectra of upscattered photons for isotropic scattering (eq. jl2| : 
dashes) andfor the case of a parallel beam of target photons (anisotropic scat- 
tering, eq. [BP: solid lines) for 7 = 10 5 and ei = 1(T 7 . In the latter case spectra 
are shown for several angles between the source and the observer, £ = (0.2) 2.2 
radians. The emissivities shown are divided by 0.002tt^. 



( is the LS angle between the momenta of the electron and incoming photon (£ = for head-on collisions, see Fig. [T]), and 

e2max = 4e 1 7 2 /(l+4e l7 ) (10) 

is the maximal energy of the upscattered photons. 

Eq. (||) gives the spectrum of the upscattered photons per single electron and single target photon, and is the basic formula in our 
derivation. For an isotropic distribution of photons this can be compared with the classical Jones' (1968) result. Let the z-axis be 
(anti-) parallel to the electron momentum so that ( = 6 and c/f2 7 = 2nd cos £; Q 1 = for the isotropic photon distribution. 
Changing the integration variable from £ to e\ gives 



d*(7,ei) 

d€2 



2efr(7- 



The result, after some rearrangement, exactly coincides with the well-known formula (lones 1968) 



dRu 



d€2 



: (7, £ i) : 



27nj 
ei7 2 



2q'\nq' + (l+2q')(l-q') + 



'\2 



1 (4ei 79 '> 



2(l+4e l79 ') 



T M-q') 



(11) 



(12) 



where o' = £2/[4ei7 2 (l -£2/7)] and 1 /47 2 < q' < 1. 

Fig. g compares the spectra of upscattered photons from a parallel photon beam (using eq. [||]) with the isotropic case (eq. p.1 
The spectra are shown versus £2/e2max for several angles ( (in the anisotropic case). The other parameters are: 7 = 10 5 and e\ = 10~ 7 . 
The spectra are different in different directions, but all are peaked at the energy which corresponds to highest possible energy of the 
upscattered photon. The difference between the two cases is maximal when f w 0. 

The interpretation of the effect is straightforward and becomes obvious for the geometry with the source and the observer located 
at the same point (£ = 0). In this position the observer sees mostly those photons which suffered head on collision and have been 
scattered back, and thus have the highest possible energy £2 = £2max- In the case of 'isotropic' scattering the number of such photons 
is negligible. Correctly averaged over the viewing angles eq. (ph will give again the same result as eq. ( 12 ). 



2.2. Integration over the solid angle fl y 

In principle, eqs. (|l]),(§), and (||) are sufficient to compute the anisotropic ICS for any photon field. However, the computation can 
be facilitated using analytical limits on the angles over which to integrate. 

We distinguish kinematic and geometrical restrictions on the angles involved. These are two independent conditions which should 
be satisfied simultaneously. The kinematics gives the range of angles £ < (0 for which a photon of energy ei can be upscattered 
into £2- For the case of a disk source the geometry defines a range of angles so that only the disk photons are considered. For both 
conditions we derive analytical formulas for the integration limits. 



4 ANISOTROPIC INVERSE COMPTON SCATTERING IN THE GALAXY 

2.2.1. Kinematic limits on Q 

Using eq. (^) gives the lower limit on cos£ 

cosC>cosCo = ^ z r-1. (13) 

2ei7(7-e 2 ) 

In terms of polar and azimuthal angles of the electron (k, tp) and incoming photon (0, cf>) (see Fig. |l]) 

cosC = - cos a« cos # + sin ac sin #cos(0 -?/>). (14) 

Eqs. ( pj| ) and ( |l4| ) place restrictions on possible combinations of the photon angles (9, </>): 

7T — K— Co < $ < I" - k + Coi 
^-arccosT < < ■i/' + arccosT, (15) 

where 

y _ cos Co + COS AC COS g 

sin aj sin 9 



In the case I T I > 1 , 



ip-ir < tfi < ip + ir, forT<-l; 

forT>+l. 



(17) 



2.2.2. Geometrical limits for a disk source 

In the case of a thin disk we can derive analytically the integration limits (eq. [^(|), which come from the position of the electron 
above the disk plane. We consider a coordinate frame centered on the electron projection on the (jc,y)-plane of the disk and with z the 
altitude of the electron above the disk (Fig. |lj). Then the coordinates of any point {x,y) in the disk plane can be expressed in terms of 
angles (9, </>): x = ztan^cos^ and y = ztan6*sin</>, where the angles 9 and cj) are essentially the same as in eqs. (Q)-(||). From the disk 
equation (x+p) 2 +y 2 < R 2 d one can derive 

(ztan6>sin^) 2 <^-(ztan6»cos0 + p) 2 , (18) 

where (x = -p,y = 0,z = 0) are the coordinates of the disk center (p > 0), and Rj is the disk radius. From this it immediately follows 
that 

c\ < cos9 < C2\ 



(19) 



Cl=z/Vz 2 + (p + Rd) 2, C2= ( z/Vz 2 + (p-R d )\ forp>R d ; 

I i , tor p K/j\ 



where 
and 

In the case |t| > 1, 



The actual integration solid angle in eq. (|8|) is, therefore, the intersection of the solid angles given by both conditions, eqs. (|15[)-(|17|) 
and eqs. 

For convenience, the angles n and can be expressed in terms of the distance of the observer from the center of the disk, R s , and 
the angle £ between R s and p relative to the disk center: 

n = ir- arctan(a / z) , 

( p-7?,cos^\ 
7r + arccos I , for < 4 < 7r, 

i>={ > n tl ( 23 ) 
7T- arccos , for -ir < £- < 0, 



(20) 

!>i=arccos-r; 02 = 2n- arccos r, 

T = (R 2 i-p 2 -z 2 t<m 2 9)/(2zptan9). (21) 



0<^<27r, forr>+l; 

(be 0, forr<-l. { ' 



where a = y (R s - pcos £) 2 + (p sin £) 2 
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2.2.3. Angular distribution of the background photons 

The angular distribution of the background photons depends essentially on the disk properties. We distinguish three extreme cases: 
transparent disk (a collection of points sources, each of which emits isotropically: Galactic plane), emitting surface (which emits 
oc cos 9: stellar surface, accretion disk), and a hypothetical intermediate case of an 'isotropic' disk. (The hypothetical 'isotropic' disk 
is considered here, since it provides an isotropic, within the solid angle fi~, distribution of background photons.) One can easily 
show that in these cases the angular distribution of photons at some particular pointQ (x = Q,y = 0,z) above the disk can be expressed 
as follows 

. / , \ -i f cos -1 9, for the transparent disk; 

2 7 (fi 7 ) = -^(0,0) ( -— 1 ) x I 1, for the 'isotropic' disk; (24) 

de 1 \de 1 J y cos9, for the emitting surface, 

where the function dq 7 {9,(j))/dej is the differential emissivity (stars and dust) of the disk (i.e., number of photons emitted per 
unit time per unit area per unit solid angle per unit energy interval) as seen from the point (0,0, z). For the case of a homogeneous 
distribution dq^/dc- = const. If the distribution of the emissivity in the plane is symmetrical about the disk axis then dq 1 {9,4>)/de 1 = 
dq 7 (r[9, (/)])/ de 7 is a function of the distance from the axis expressed in terms of the angles, r(9,<j)) = (p 2 +z 2 tan 2 6*+2pztan6'cos</>) 1 / 2 . 
The differential number density of background photons at any point (0,0, z) is given by 

dHl - ^dcose / d^A(cos0,<£>x { 1 }. (25) 



de 7 J Ci J fa de 7 j cos g 

where c\ ,2,^1.2 are defined in eqs. (^o|)-(^2|). For the simplest case dq 1 /de 1 = C = const, at a point slightly above the disk center 
(4>i = 0,4>2 = 2tt,c\ — > +0,C2 = 1) one can obtain 

dn { lo S s 1 S=C1 { - 27rC1 °g c l 1 

p=C(^2-0i)x { s ) ={ 2ttC \, (26) 

where the infinity in the top row follows from using the thin disk appr oxim ation and does not create any singularity in the formulas 
since it appears only for points at the disk plane, z = (see also Section 2.4 ). 



2.3. Scattering off a spectral distribution of electrons and source photons 

The final 7-ray emissivity spectrum (photons cm" 3 s -1 sr" 1 MeV" 1 ) can be obtained after integration over the energy spectra of 
electrons and the background radiation 

dF n e n 7 c f f 2 2 x <iR(7,ei) 



. - ,,4, djefYM^feh) V (27) 

dt2 Aitm e c L J J de-2 

where (n 7 e 2 / 7 ), (n c 7 2 / e ) are the differential number densities of photons and electrons, respectively. The resulting spectrum depends 
strongly on 62, / 7 , e , and spatial coordinates. However, we find that for the power-law electron spectrum, at least, the ratio 

^Ptt (28) 

dF iso /d€2 

is much less sensitive and depends mostly on the spatial variables. The value of dFj SO /de2 is calculated in the same way as eq. (p7|), 
but instead of dRjdt2 the isotropic function is used dRj SO /de2 (eq- JT^])- Therefore, for particular geometry and photon and electron 
spectra, the value 3? gives essentially the geometrical factor. 

2.4. Calculations and Discussion 

The effect described has a clear geometrical interpretation. Jones' (1968) formula assumes an isotropic distribution of photons 
while in our case the photons are concentrated in a solid angle £7 7 which is defined by the geometry. Therefore, 47r/f2 7 gives a 
rough estimate of the enhancement when the effect of anisotropy is taken into account. This is particularly true when one considers 
the maximum energy of the upscattered photons £2 = £2max(7,£i) which can be reached for given e\ and 7 (eq. [|I(i|]). At energies 
£2 < £2 max, this effect is smaller since the angular distribution of the background photons becomes less important (larger scattering 
angles in the ERS become allowed). At very small energies t2 *C £2max> the effect disappears completely (anisotropic = isotropic) 
since this corresponds to the case when the background photons are scattered at large angles, and their actual angular distribution is 
of no importance. 

Fig. H shows the spectral emissivity ratio J— /^f of anisotropic to isotropic scattering (eqs. [||, [|l2||), vs. z for electrons at the disk 
axis (p = 0) for an observer in the center of the emitting disk (left), and at the solar position R s = 8.5 kpc (right). We show three cases: 

2 Note, we consider the coordinate frame where the (x,y)-p\ane coincides with the disk plane and the z-axis passes through the position of the electron, not the disk 
center! 
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FIG. 3. — The ratio -j^/^j^- (eqs. Jl^]) vs. z for a small volume of isotropically distributed electrons at the disk axis for an observer at the center of the 
Galaxy (left), and at the solar position R s = 8.5 kpc (right). Dotted lines show the ratio for the 'isotropic' disk, solid lines for the transparent disk (the Galaxy), and 
dashes for the emitting surface case. Other parameters: Rj = 15 kpc, e\ = 2 X 1(T 7 , 7 = 10 5 . 



transparent disk (the Galaxy), emitting surface, and 'isotropic' disk. The other parameters are: Rj = 15 kpc, t\ = 2 x 10~ 7 , 7 = 10 5 . 
To illustrate the energy dependence of the ratio, the curves are shown for three different photon energies: 62 ~ £2max, £2max/l-2, 
and 62 m ax/ 1-5. In the simplest geometry (left panel) when the scattering electrons and the observer are both on the disk axis, the 
'isotropic' disk case gives exactly the same enhancement as obtained from the rough estimate 47r/f2 7 . This reflects simply the larger 
number of background photons (by a factor of 47r/51 7 ) per solid angle compared to the isotropic case. Other cases, the transparent 
disk and emitting surface, give somewhat different results because of the more complicated angular distributions of the background 
photons. The ratio increases with z until the whole disk is covered by the solid angle f2 7 , which is defined by eqs. (15)-(p^), and 
then remains almost a constant. 

The importance of the angular distribution of the background photons is clear from comparison of the behaviour of the ratio at 
small z on the left and right panels (Fig. |3|). In the case of the transparent disk, the photon distribution is ~ cos" 1 9 which has a 
minimum at 8 = 0. The solid line in the left panel (scattering electrons and the observer are both on the disk axis) thus shows the 
smallest effect. The picture changes when the scattering electrons do not lie just above the observer (right panel). In this case, z — > 
corresponds to 6 — > ir/2 where the angular distribution has an infinite peak. Therefore, in the case 62 ~ £2max, the ratio — » 00 as 
z — ► 0. This infinity, however, is not physical since it is connected with approximation of infinitely thin disk plane. The result can be 
compared with the emitting surface case (~ cos 6*; dashed line) where the opposite, but finite, effect is seen. 

Fig. @ shows 5ft (eq. [§§]) as calculated for a homogeneous transparent disk with radius Rj = 15 kpc, (.2 = 100 MeV, monoenergetic 
background photons e\ = 1 eV, an electron spectral power-law index of -3, and for z = 1,4, and 10 kpc. At small distances from the 
plane within the disk the ratio is close to unity since the emitting electrons directed to the observer are moving under small angles to 
the disk plane and the radiation field they see is effectively isotropic. As the ICS electrons become higher in altitude and/or further 
from the disk center the radiation field they see becomes more and more anisotropic, and therefore the ratio 5ft increases. 

The smaller value of the effect compared to that shown in Fig. || is the result of integration over the electron spectrum. As 
mentioned above the effect is maximal when £2 = £2max or equivalently 7 = 7 m i n (ei,e2) = \(^2 + [^ + £2/ 'ei] 1 ^ 2 )- At 7 > 7 m ; n the 
effect is smaller, but the background photons come from a larger solid angle to be upscattered to the same energy; this increases the 
number of background photons involved. The actual value of the enhancement is then naturally obtained from a balance between the 
decreasing number of electrons at higher energies and the increasing value of the effective solid angle allowed for the background 
photons. 

Fig. Q (bottom right) shows as an example a calculation made for the emitting surface (z = 1 kpc). An enhancement of the emission 
above the solar position is clearly seen. 

3. DIFFUSE GALACTIC 7-RAYS 

In the following the reader is referred to SMR99] for details of the models. 

The interstellar radiation field (ISRF) is essential for propagation and 7-ray production by CR electrons. It is made up of contri- 
butions from starlight, emission from dust, and the cosmic microwave background. Estimates of the spectral and spatial distribution 
of the ISRF rely on models of the distribution of stars, absorption, dust emission spectra and emissivities, and is therefore in itself a 
complex subject. 

Recent data from infrared surveys by the IRAS and COBE satellites have greatly improved our knowledge of both the stellar 
distributi on and the dust emission. A new estimate of the ISRF has been made based on these surveys and stellar population models 
flSMR99| ). Stellar emission dominates from 0.1 fim to 10 fim, emission from very small dust grains contributes from 10 fim to 30 
fim. Emission from dust at T ~ 20 K dominates from 30 fim to 400 fim. The 2.7 K microwave background is the main radiation 
field above 1000 /zm. 

The ISRF has a vertical extent of several kpc since the Galaxy acts as a disk-like source of radius ~15 kpc. The radial distribution 
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FIG. 4. — Anisotropic/isotropic ICS ratio !R (eq. ]28fl) as calculated for an observer at R s = 8.5 kpc from Galactic center. Other parameters are: Rj = 15 kpc, electron 
power-law index -3, ej = 1 eV and €2 = 100 MeV, for z = 1 kpc (top left), 4 kpc (top right), and 10 kpc (bottom left). The bottom right panel shows a calculation for 
the emitting surface case (z = 1 kpc). 



of the stellar component is centrally peaked because the stellar density increases exponentially inwards with a scale-length of ~2.5 
kpc until the bar is reached. The dust component is related to that of the neutral gas (HI + H2) and is therefore distributed more 
uniformly in radius than the stellar component. 

In practice we calculate the anisotropic/isotropic ratio K for any particular model of the particle propagation (halo size, electron 
spectral injection index etc.) on a spatial grid taking into account the diff erence between stellar and dust contributions to the ISRF, 
and then interpolate it when integrating over the line of sight (see |SMR99 ). 

Fig. [5] shows a Galactic latitude - longitude plot of the intensity ratio for 11.4 MeV 7-rays for two Galactic models with halo 
size zi, = 4 kpc and 10 kpc. This is obtained from the computed sky maps in the anisotropic and isotropic cases. The calculation 
has been made with a 'hard' interstellar electron spectrum (the interstellar electron spectrum is discussed below). It is seen that the 
enhancement due to the anisotropic ICS can be as high as a factor ~1.4 for the pole direction in models with a large halo, Zh J> 10 kpc. 
The maximal enhacement occurs at intermediate Galactic latitudes in the outer Galaxy because there head-on collisions dominate. 
Fig. H shows the intensity ratio vs. 7-ray energy for several directions as seen from the solar position for a halo size Zh = 4 kpc and 
10 kpc. 

Fig. f7| shows the spectra of the inner Galaxy {\b\ < 5°, 330° < I < 30°) and for high Galactic latitudes ( \b\ > 70°, all longitudes) as 
calculated for the HELH model. The HELH model ('hard electron spectrum and a large halo', SMR99( ) is chosen to fit high-energy 7- 
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FIG. 5. — Latitude - longitude plot of the anisotropic/isotropic intensity ratio for 11.4 MeV 7-rays (ratio of the two sky maps). Halo size z^ = 4 kpc (left) and 10 
kpc (right). 




FIG. 6. — The intensity ratio vs. 7-ray energy for some direction as seen from the solar position. The corresponding Galactic coordinates (Lb) are shown near the 
right scale. Halo size Zh = 4 kpc (left) and 10 kpc (right). 



rays using hard electron spectrum^ (the injection spectral index is taken as -1.8, with reacceleration) and a broken nucleon spectrum 
(injection spectral index is -1 .8/-2.5 with break at 20 GeV/nucleon, with reacceleration). Such a nucleon spectrum satisfies the limits 
imposed by antiprotons and po sitrons ( MSR98| , |SMR99 ). The halo size taken, Zh = 10 kpc, is within the limits (4-12 kpc) obtained 
from our ll) Be studies (SM98). Following Pohl & Esposito (1998), the consistency with the locally measured electron spectrum 
above 10 GeV is not required, since the rapid energy losses cause a clumpy distribution so that this is not necessarily representative 
of the interstellar average. For this case, the interstellar electron spectrum deviates strongly from that locally measured as illustrated 
in Fig. ||. Because of the increased ICS contribution at high energies, the predicted 7-ray spectrum can reproduce the overall intensity 
of the inner Galaxy from 30 MeV to 10 GeV. If the hypothesis of a hard electron injection spectrum is correct in explaining the large- 
scale Galactic emission, then the same model for the electrons is applicable to the halo where large-scale propagation is relevant. We 
consider a wide high-latitude region 70° < b < 90° which encloses a large volume of halo emission and is therefore unaffected by 
short-scale fluctuations. 

3 In fact, the size of the anisotropic effect will increase for a softer electron spectrum (see Section |2.4| ). 
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FIG. 7. — Gamma-ray energy spectrum of the inner Galaxy (left panel) and from high Galactic latitudes (\b\ > 70°, all longitudes) (right panel) as compared with 
calculations in the HELH model. Separate comp onents show the, contrib ution of hrem sstrahlimg. 7r°-riec ay. and ICS with (solid line^ and without (do tted line) the 
anisotropic effect. Data in the left panel: EGRET ( Strong & Mattox 1 gggl COMPTF.T ( [strong etaL I9qq|> and Q SSE ( Kinzer Pnrcell Kurfess IQQc VData in the 
right panel: COMPTEL hi gh -latimde total intensity (open boxes: BloeTnen et al. 1999 diamonds: [Kappadath 1998 crosses: Weidenspointner et al. 1999] ), and total 
from EGRET measurements ( 5MR99 1. 



cv 
I 



i 

is 



10 



Interstellar electron spectrum 



30-5pMeV 1.0(1(^80. Q / /181.0(1< | 359.0 




10" 10^ 10" 10^ 

Electron momentum, MeV/c 



1 . 4 - galprop 



S 0.8 



0.0 




G a I a ctic Latitude 



FIG. 8. — Electron spectra at R s = 8.5 kpc in the plane, for 'conventional' 
(solid line) and hard electron spectrum, (dashes) models. Data points: direct 
measurements, see references in MS98E . 



FIG. 9. — The latitude 30-50 MeV 7-ray profile of the Galaxy as calcu- 
lated in the HELH model (heavy line; convolved) compared to EGRET data 
(dots; point sources removed). Separate components show the contribution 
of bremsstrahlung (thin histogram), 7r°-decay (thick histogram), and ICS with 
(dashes) and without (dash-dots) anisotropic effect; horizontal line: isotropic 
background. 



For comparison the dotted lines in Fig. [7] show the spectra of ICS 7-rays as calculated using the isotropic Jones' formula (eq. Jl2|]). 
The model with anisotropic ICS component gives somewhat larger fluxes. The effect is larger at high Galactic latitudes, but also 
significant for the inner Galaxy. For the case of the inner Galaxy the result is not sensitive to the halo size and for smaller halos the 
plot would be almost the same. At high latitudes, where the ICS is the only important contributor into the Galactic diffuse emission, 
it provides up to a factor ~ 1 .4 larger flux compared to isotropic calculations. 

Fig. U shows the 7-ray latitude profile of the Galaxy for 30-50 MeV as calculated in the HELH model, again with and with- 
out the anisotropic effect, compared to the EGRET data. For this comparison the 7-ray sky maps computed in our model have 
been convolved with the EGRET point-spread function; point sources contribution has been removed from the EGRET data (see 
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SMR99 for details). At high latitudes the model with anisotropic ICS component again provides significantly larger fluxes com- 
pared to isotropic calculations. This is particularly important for estimates of the extragalactic background since Galactic ICS and 
extragalactic contributions may be comparable. 

4. CONCLUSION 

The ICS is a major contributor in the diffuse Galactic emission above ~ 1 MeV. Since accurate observational data from EGRET is 
now available it is desirable to compute the ICS including all important effects. We have shown that the anisotropy of the interstellar 
radiation field has a significant effect on the intensity and angular distribution of 7-ray radiation, and can increase the high-latitude 
Galactic 7-ray flux up to 40%. This effect should be taken into account when calculating the Galactic emission for extragalactic 
background estimates. 

The der ivations presented here provide the basis for the treatment of ICS in our companion study of Galactic diffuse continuum 



emission (SMR99). 
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